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ON THE MICROLOCAL ANALYSIS OF THE GEODESIC X-RAY 
TRANSFORM WITH CONJUGATE POINTS 

SEAN HOLMAN, GUNTHER UHLMANN 


Abstract. We study the microlocal properties of the geodesic X-ray transform A on a 
manifold with boundary allowing the presence of conjugate points. Assuming that there 
are no self-intersecting geodesics and all conjugate pairs are nonsingular we show that the 
normal operator JV = A ( o X can be decomposed as the sum of a pseudodifferential opera¬ 
tor of order —1 and a sum of Fourier integral operators. We also apply this decomposition 
to prove inversion of X is only mildly ill-posed in dimension three or higher. 


1. Introduction 

Our object of study in this paper will be the geodesic X-ray transform defined on a 
Riemannian manifold (M. g ) with boundary. Loosely speaking, the geodesic X-ray trans¬ 
form, which we denote in the present work by A, is the operator which takes a function 
defined on M to its integrals along all geodesics of (. M,g ). The Euclidean version of this 
problem provides the mathematical basis for X-ray computerised tomography, and has a 
long history going back at least to the well known work of Radon [22] . Slightly before 
Radon, Funk also found an inversion formula for the case of symmetric functions on the 
two-sphere [5]. Study of the general non-Euclidean case began in earnest, to the author’s 
knowledge, with Mukhometov m in relation to the boundary rigidity problem which is 
the problem of determining a Riemannian metric from knowledge of its distance function 
restricted to the boundary. The boundary rigidity problem arises in seismology when we 
take the so-called “travel-time metric” for seismic waves, and consider that we can measure 
the amount of time it takes for a seismic wave to travel from a source to a receiver location. 
Indeed, the geodesic X-ray transform of a tensor field arises as a linearization of travel time 
tomography [28] . Other applications of the geodesic X-ray transform include ultrasound 
transmission tomography m and optical tomography with a variable index of refraction 
(12j. It has also been shown that injectivity of the geodesic X-ray transform can imply 
identifiability results for the anisotropic Calderon’s problem in some special cases [3]. 

Most existing results on the geodesic X-ray transform for manifolds with boundary 
concern simple manifolds. A simple manifold (M,g) is one for which the exponential map 
centered at any point is a diffeomorphism onto M, and dM is strictly convex. Using the 
energy integral method originally introduced by Mukhometov [13], it can be shown that 
the geodesic X-ray transform is injective, and stability estimates for the inversion can be 
obtained, under hypotheses on the sectional curvature which imply that (. M,g ) is simple 
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(see [231). ^ has also been known for some time that the normal operator 

AT = X t o X 

is an elliptic pseudo differential operator of order —1 when ( M,g ) is a simple manifold 
[9]. The transpose X t must be defined by placing an appropriate measure on the space of 
geodesics. An explicit Fredholm type formula which can be inverted by a Neumann series 
in the case of a simple manifold was found in the two dimensional case in [20]. The geodesic 
X-ray transform acting on tensor fields has also been studied extensively almost completely 
in the case of simple manifolds. See [19j and the references therein for a survey of recent 
progress, and also a listing of current open problems. See [18] for the most general existing 
results on tensor tomography for simple surfaces. 

There are few results for non-simple manifolds and most of those that exist are relatively 
recent. By applying analytic microlocal analysis injectivity of X can be proven in some 
non-simple cases [11126]. One of the most general results to date for dimension three or 
higher was established using the scattering calculus and a layer stripping argument ]29], and 
applies in some non-simple cases. It is difficult however to connect the hypotheses in [29] 
to other more geometric assumptions on ( M,g ). Recent work using tools from dynamical 
systems also shows that X, possibly acting on tensors, is injective for manifolds with 
hyperbolic trapped sets, although not including any conjugate points [6]. Similar results 
for the transform of a connection are shown in [7]. In the two dimensional case injectivity of 
the geodesic X-ray transform was established for some non-simple cases including conjugate 
points in [24] . although this injectivity should be held in contrast to other recent work 
showing that when there are conjugate points present in two dimensions it is not possible 
to establish a stability estimate for inversion of the geodesic X-ray transform between any 
Sobolev spaces ns]. Some results for the tensor problem on a non-simple manifold have 
also been found in [2]. 

In general, the normal operator Af is not a pseudo differential operator when (M, g ) is not 
simple. The authors of m show that in the case of fold caustics an appropriately localised 
version of the normal operator is the sum of a pseudodifferential operator and a Fourier 
integral operator. This result is very much in the same spirit as the results of the current 
paper, although here we lessen the restriction to fold caustics. In m the restriction to fold 
caustics is also removed although only in the case of two dimensions. In fact the method of 
m is similar to the current paper, but here we have analysed the geometry of conjugate 
points in more detail in order to reach a more general conclusion. 

This paper contains two main results. The first is Theorem [4] which shows that when 
there are no singular conjugate points the normal operator A f can be decomposed as an 
elliptic pseudodifferential operator of order —1 plus a sum of Fourier integral operators 
(FIOs). Each FIO corresponds with conjugate pairs of a given order, and the order of the 
FIO depends on the dimension of M and the order of the conjugate points. The canonical 
relation of each FIO depends on the geometry of the conjugate pairs. The theorem is 
actually a bit more general than this, allowing for a weight to be included in X. The 
second result is Theorem [5] which concerns the use of the decomposition in Theorem [4] to 
obtain stability estimates for inversion of X. While results similar to Theorem [5] have been 
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shown before we are not aware of this exact form which we have included because it matches 
with hypotheses in previous literature on the the convergence of Tikhonov regularisation. 
This is described in more detail in a remark after the statement of the theorem. 

2. Preliminaries 

In this section we introduce the notation which will be used throughout the paper, a 
few important definitions, and some preliminary lemmas we will use later. Throughout the 
paper we make the following assumption on (. M,g ). 

Assumption 1. ( M,g ) is an n dimensional compact, Riemannian manifold with smooth 
strictly convex boundary and with n > 2. Further assume that ( M,g ) does not contain any 
self-intersecting geodesics. 

The requirement that ( M,g ) does not contain self-intersecting geodesics is stronger than 
the requirement that ( M,g ) be non-trapping (i.e. not contain any closed geodesics). The 
prohibition against self-intersecting geodesics arises for technical reasons in the application 
of the calculus of FIOs, and may be able to be lessened to simply non-trapping via an 
additional layer of microlocalization, but we have not done this. 

The interior of M will be M vnt = M \ dM. The unit sphere bundle for M (resp. M mt ) 
will be SM (resp. SM mt ). We will use it for the projection mapping ir : SM — > M, and 
ism for the inclusion map isM '■ SM —>• TM. Further, we will use the same notation for 
7i \sm int and isM\sM int (this should not cause any confusion). We will also have occasion 
to require other bundle projection mappings, and these will always be denoted by n with 
different subscripts indicating the bundle. For example, for the projection from T*SM mt 
(this is the cotangent bundle of SM int ) to SM int we will write it T * SM int. 

In referring to points of vector bundles a timeless question is whether the notation should 
include the base point. For example, when indicating points in the tangent bundle TM 
should one use ( x,v ) £ TM or just v £ TM7 We do not fully commit ourselves to either 
choice, but stipulate that when we write ( x,v ) £ TM this is actually a short hand for 
v £ TM and x = n The same holds for points in other vector bundles. 

At the boundary dM, the outward pointing unit normal vector will be u + , and the set 
of inward or outward pointing unit vectors will be 

d ± SM = {(x,u) £TM : x£ dM, ±{v,u + ) > 0}. 

Note that d-SM is the set of inward pointing vectors, and d+SM is the set of outward 
pointing vectors. 

For v £ TM we write 'y v for the maximally extended geodesic with initial data A f v (0) = v. 
The reader should specifically note here that we consider 7„(0) to be a point in TM, and 
so, loosely speaking, this initial condition is including both the initial position and tangent 
vector of the geodesic. 

We define two functions t± initially on SM mt by the formulas 

7r(/r±(T)) e dM, 
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t+(v) > 0, and r_ (v) < 0. These give the backward and forward “time to the boundary.” 
Since the boundary is strictly convex t± £ C°°(SM int ), although t± are not smooth up to 
the boundary [23]. Nonetheless r± extend to smooth functions on SM \ TdM, and thus 
we have t± £ C°°(SM \ TdM). 

We will make much use of the smooth mappings F : SM mt —> d-SM defined by 

F(v) = 7 V (r-(v)), 

and T : V = {( v , s ) £ SM int xl : r_(u) < s < t + (u)} -a- SM mt defined by 

( 1 ) d>(v,s) =%(s). 

The mapping \k is the geodesic flow in the unit sphere bundle, and F takes v £ SM int to 
the unique w £ d-SM such that the tangent vector to the geodesic with initial data w 
equals v at some point along the geodesic. We will also write \k for the geodesic flow on 
TM. 

On M we have the Riemannian density which we write as |du 9 |. Further, g induces 
Liouville densities on SM and dSM which we will write as |dSM| and |d<9SM| respectively. 
With these notations the Santalo formula (see [23]) is 

r r r r +( w ) 

(2) / h(y) |d5M(u)| = — / / h{^/ w {s)) (w^v + ) ds \ddSM(w)\ 

J SM J d-SM Jo 

which holds for h £ C{SM). In view of this formula we use |d/Li| = —(rc,u + ) IdcIS'Ml as 
the density on d-SM. 

1/2 

When dealing with FIOs we need to use half densities, and so we will write 0^ for 
the half density bundle on M. The same subscript notation will be used for half densities 
on other manifolds such as ihj SM for the half densities on d-SM. We will also consider 

1 /2 

function spaces of sections of half-densities writing for example C£° (^ Afrnt ) for the space of 

sections of fwhich are smooth with compact support contained inside of M int . We now 
define the main object of our study, the geodesic X-ray transform, acting on half densities. 

Definition 1. (Geodesic X-ray transform on half densities) The geodesic X-ray 
transform on ( M,g ) is defined for f £ Cf°{Ttas a half density on d-SM by the 
formula 

X[f]{v) = | d /|i/ 2 (7^( g )) |dM(^)| 1/2 - 

It can be shown that X : Cf° (Ilj/L*) —X Cffi(Tl l Q 2 SM ) continuously. In fact X is an FIO and 

can be extended to distributions £'(Q^‘l at ) and various other spaces (including L 2 (Q^ 2 . nt )). 

Central tools for our analysis will be the push-forward and pull-back of half densities 
which we now define. 

Definition 2. (Push-forward and pull-back of half densities) Let X and Y be two 

manifolds of dimensions nx and uy, and let G : X —x Y be a smooth submersion. Also, 
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suppose we are given half densities |d^x| 1//2 and |d/iy| 1//2 on X and Y. The push-forward 
G* : Cf°(X,n 1 / 2 ) -a G c °°(^y /2 ) and pull-back G* : G c °°(Qy /2 ) -A G°°(Q^ /2 ) are defined by 
the requirement that 

f y h(y) G*[f](y) = J x | d ^|i/ 2 ( G, ( x )) \ d Tx\ 1/2 f(x) = lx G*[h](x) f(x) 

for all f € G c °°(^ /2 ) and h € C c °°(I^/ 2 ). 

It can be seen directly from this definition that the pull-back and push-forward are adjoints 
of one another. Intuitively the pull-back is the precomposition with G , while the push- 
forward is an integration over the level sets of G. 

Indeed, let us consider the push-forward and pull-back by the map F : SM int — > d-SM 
introduced above using the half-densities corresponding to the densities \dSM\ and |d//| 
also introduced above. Using Santalo’s formula Q we find that 

[ Hw) F* [/] (w) = [ -L—(F(v)) \dSM(v)\V 2 f(v) 

Jd-SM JSM |d/i| 

= Ssm xL (m) ]dsm^ (v) lASM(v)l 

= Lm (T'” wm» iUe)) ds ) wy ( w) |d ' ,wl 

= Lm {L +< ’ wk^ iU ‘ }) dS ) Id " (, " )|1/2 HW) 

From this we have explicit formulae 

( rT +( w ) r \ 

Jo \dSM\V 2 LL S )) ds j \dp,(w)\ 1/2 

and 

(4) F*[h](v) = -*(F(v)) |dS , M(u)| 1//2 . 

Similar analysis, not requiring the Santalo formula, may be applied to it : SM mt —» M to 
find that 

n,lfl{x) = {Is,M \dSM\^ (V) |dS * MWI ) |d ” s(l)|1/2 

and 

^*[h]{y) = \M v)) \dSM(v)\ 1/2 . 

\Av g \F^ 

Using these formulae and the observation above that the adjoint of a pull-back is a push- 
forward, we find that the geodesic X-ray transform and its adjoint may be written as 

X = F* o tt* and X t = 7r* o F*. 
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The normal operator is thus 

Af = 7T* o F* o F* o 7r*. 

This will be the jumping off point for our analysis of A f as an FIO. 

We actually generalise slightly by adding a weight. Given £ C°°(SM ) we will write 
(f) m for the operator which multiplies by </>, and define 

X$ = F* o cj) m o n* 

This is the weighted geodesic X-ray transform, and since it requires little extra effort to 
prove our results for X^ we will do so. 

We will now show that it is generally true that the push-forward and pull-back by a 
submersion are both FIOs. This is certainly a known result going back to at least [8], but 
we provide it here for completeness, and so that we have it in precisely the form required 
for the rest of the current paper. 


Lemma 1. Suppose we are in the setting of definition [H Then the pidl-back G* and the 
push-forward G* are both FIOs of order (ny — nx )/4 with non-vanishing principal symbol. 
The canonical relation of G* is 

C G . = {(£, DGH 0 : x £ X, £ £ T* (x) Y \ {0}} 

while the canonical relation of Cq* is 

G g . = {(DG|U,0 : x£X, ^T* (x) Y\{0}}. 

i—i 1/2 

Proof. Let Kq* be the Schwarz kernel of G*. From definition [2], for any / £ C£° (fi) and 
h £ G£°(fiy 2 ) we have 

f h f 

(K XsM ,h ®f) = j x |d ^ [1/2 (G(x)) |d ^ |1/2 (^) \dfix(x)\ 

Suppose now that the support of / and the support of h are contained in the domains of 
single coordinate charts <f> : X —> M nx and <j> : Y —>• M ny respectively, that 

G = fioGofi- 1 : M nx M ny 


is the representation of G in these coordinates, that 

J(x) |da;| = </>*( \dnx\) 

is the representation of the relevant density in the coordinates on X, and that 


/(*) = 


/ 


1 (x)), h(y) = 


\dfi X \ 1/2 

Then using the Fourier inversion formula 


h 




|d^y| 1/2 


( 5 ) 


( K x SM i h ® f) = / HG(x)) }{x) J(x) dx 

JR n x 


e l£-(G{x) dy 


(27r) n v 


lj y x*" y xM; x 
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Untwining the definitions we can see that this is a local representation of the Schwartz 
kernel of an FIO of order (ny — n x )/4. From the phase of the local representation ([5|) we 
can further see that the claimed canonical relations are correct. Note that since DG\ V is 
a surjective map at every v, DG |* is injective and so in fact Cq* and Cq* are canonical 
relations. Finally, from the local representation ([5j) we see that the principal symbol is 
equal to J(x) multiplied by an appropriate non-vanishing half density, and so the principal 
symbol also does not vanish. □ 

Since the dimensions of M, d-SM and SM are n, 2 n — 2, and 2n — 1 respectively, Lemma 
|T] shows in particular that 

F* £l~ l/4 (d.SMxSM, C' K] Oy_ 2 5MxSM ), F* erV^SMxd-SM, C' F rt'* SMxSM ), 

7T* € X (1 " n)/4 (M x SM, Cnf xSM ), and € 1^\SM x M, ^ xM ). 

The bulk of the remaining analysis on the normal operator J\f will be to show how the 
clean composition calculus may be applied to analyse the compositions of these operators. 
First however, we review some more preliminary material. 

The canonical symplectic form on T*M will be denoted uj, and the corresponding one 
on TM induced by g denoted by oj g . Specifically, if \> g : TM int —> T*M int is the musical 
isomorphism given by g then uo g is defined by 

u g (X,Y) =u(D\> g X,D\> g Y) 

for any X and Y in T(TM mt ). This symplectic form then induces an isomorphism 
from T(TM int ) to T* (TM mt ) at each point defined by 

[b Ug {X)\{Y)=u g {X,Y) := i x uj g (Y) 

for all X and Y in T(TM int ) (here \ x is interior multiplication by X which is defined by 
the previous formula). We denote the inverse of this map in the usual way by (L g • 

The geodesic flow on TM is given by the Hamiltonian H{v) = l/2(||u||^ — 1) € C°°(TM). 
Thus we have for v £ TM int the following invariant formula for 7„(0) € T(TM int ) 

%(0) = X%(v) 

where X F 9 ( v ) is the Hamiltonian vector field given by H with respect to the symplectic 
form uj g evaluated at v. Since SM is a level surface for H in fact X F 9 (v) £ TSM int and 
we see that 7„(0) defines a smooth vector field on SM mt . We record these observations as 
a lemma which will be useful later. 

Lemma 2. The smooth vector field X F 9 on TM restricts to a smooth vector field on SM int 
and for each v £ SM int 


U0) = x u h 9 (v). 
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3. Characterisation of A f 

In this section we reach the main result of the paper, which is a characterisation of AQ, 
as a sum of Fourier integral operators. To do this we use the decompostion described in 
the previous section 

(6) A/0 = 7T* o (j) m o F* o F* o </> m o 7T*. 

Our first result in this direction is the following theorem. 

Theorem 1. Define 

Nsm = F* o F*. 

Then N S m G X -‘/ 2 [sM int x SM int , C' Msm ; where 

CMsm = {{DF^DFlU) : v, v £ SM int , F(v) = F(v), £ £ T* F{v) d.SAl} . 
Furthermore, the principal symbol of Ms m does not vanish. 

Remark 1. It is worthwhile to note that, while we are using the abstract FIO composition 
calculus to analyze the operators in this theorem, in fact it is not hard at all to find an 
explicit formula for Nsm ■ Indeed, using ([3|) and flU) we can see that 

N su [h]{v) = d* IdSA/l'T 

Proof. The proof will be an application of the clean composition calculus for FIOs. There 
is an adjustment required since F* is not properly supported. However, using the fact that 

F* : t C^(n^ SMint ) 

continuously and F* : ) -A- C°° there is no issue defining the compo¬ 

sition of the two operators, and we may localize and reduce to the case of a composition of 
two properly supported FIOs. Therefore the proof is reduced to analysis of the canonical 
relations. 

Using Lemma [1] it is straightforward to see that 

CjVsm = Cf* o Cf, , 

where Cjy SM is defined in the statement of Theorem HI and so if we can show that this 
composition is clean the proof will be almost complete. To show the composition is clean 
we must show: 

(i) The intersection 

(7) C = (C F * x C Ftt ) n (: T*SM int x A (T*d-SM) x T*SM int ) 

is clean in the sense that C is an embedded submanifold, and at every point c £ C 
the tangent space T c C is equal to the intersection of the tangent spaces of the two 
manifolds being intersected. 

(ii) The projection map 7rc : C —>■ T*SM int x T*SM int is proper. 
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(iii) For every (j> G T*SM int x T*SM int , 7 r c , 1 (</>) is connected. 

To begin we note that since DF\ V is surjective for every v, DF\^ is injective for every 
v. This implies that ttq is injective, and so point (iii) is immediate and assuming the 
other points are proven the excess is zero. Since the principal symbols of each of the two 
operators is not vanishing this will also show that the principal symbol of the composition 
is not vanishing, and so in fact all that remains is to demonstrate points (i) and (ii) hold. 

To begin proving point (i) we first observe that 7„(0) is in the kernel of DF\ V for all 
v € SM mt , and since DF\ V is surjective by a dimension count the full kernel of DF\ V is 
the span of 7 „( 0 ). By Lemma [2] the range of DF |* is thus 

{6 € T*SM int : 9 {X%(v)) = 0} . 

Letting v vary, and using the nontrapping assumption on ( M,g ), we find that the range of 
the map 

G : {(£, s) € T*d.SM xK : 0 < s < r+(£)} 3 (£, s) ^ DF\\ , s) £ € T*SM int 

is precisely the set 

( 8 ) {9 € T*SM int : 9 t t * S m(9))) = 0} . 

Again by Lemma [2] the map 

T*SM int 3 9^9 t t *sm(B))) 

is a smooth submersion and so the set ([5]). as the zero level set of a smooth submersion, is 
a closed embedded submanifold of T*SM mt . Further, the map G is a diffeomorphism onto 
that set as can be shown using the inverse function theorem. We will write the inverse of 
G composed with projections onto the £ and s components as G^ 1 and G 7 1 respectively 
Now note that C can be parametrized by a mapping 


P c : V Pc = {(£, s, s) € T*d-SM xlxl : 0<s,s< r + ( 7 r T . 5M (£))} -»■ C 


defined by 


Pc(Z,8,S) = (DF\ t % 


T*SM 


«)( s ) 


DF\* 


T-7T. 


T* SM 


d)G) 


0 


We initially have directly from its construction that Pc is an injective immersion into 
Range(G') x T*d-SM x T*d-SM x Range(G), and to show that it is a smooth embedding 


it only remains to show that the inverse is continuous. To see the inverse P ( , 1 is continuous 
we note that it is the restriction to C of the map 

~%ext ,—1 


p, 


c 


: Range(G) x T*d-SM x T*d-SM x Range(G) -A V Pc 


defined by 

P^\A 1 ,A 2 ,A 3 ,A^) = 

Since G is a diffeomorphism, in particular i s continuous, and so is Pq 1 . Thus Pc 

is a smooth embedding and C is a 4n — 2 dimensional embedded submanifold. 
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To check that the intersection is clean note that, similar to C, Cf* x Cf, can be 
parametrized by a related mapping 

Pc Ft xCp* = v Pc F * xCf , = {(€, l s, s ) G T*d-SM x T*d-SM x K x R : 

o < s < T + (-Kt*Sm(0)i 0 < S < T + (lTT*SM(i))} 

-> T*SM int x T*d_SM x T*d_SM x T*SM int 

defined by 


Pc F * xC F , (C) £> = ( ^1^ ,«(«)£> -(*)£ 

y , 7 r T*SM ^' k ' Sm(€) v 

Suppose that Co € C. Then any tangent vector X € T cq (Cf* x CV*) comes from a smooth 
curve 

^(a),f(a), s(a), s(a) 

such that Pc F *xC Ft (£(0),£(0),s(0),s(0)) = co via 

Pc F *xC F „ U(a),£(a),s(a),s(a)) . 




<r=0 


This tangent vector X is also in T C0 (T* SM int x A (T*d-SM) x T*SM int ) if and only if 


_c^ 

da 


£(«) = 


a=0 


da 


|(a). 


a=0 


Thus expressed in any local coordinate system £ and £ agree to first order at a = 0, and 
so in fact 


X = 


da 

_d_ 

da 

_d_ 

da 


<r=0 


Pc F *xC F , (£(a),£(a), s(a), s(a 


Pc F »xC F , (£(a),£(a),s(a),s(«)) 


<T=0 


Pc (£(a),s(a),s(a)) . 


<r=0 


Therefore X € T Co (7, and the intersection is clean. Thus point (i) is proven and all that 
remains is point (ii). 

For point (ii) we note that if K <s T*SM int x T*SM mt is a compact set, then for some 
constant e > 0 we have that 


K <s L := {(0, 9) € T*SM int x T*SM int : dist(vr o vr T , SM int ( 6),dM) > e, 

dist(7T O 7T r * S M int{0),dM ) > e, 

| 5 + ll4)<^} 
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where L is compact. Now, because the boundary of M is strictly convex there exists e > 0 
such that for ( 9 , 9) € L 

(F(7r T , SM int(9)),i' + )g < — e, (F(iT T , SM int(9)), v + ) g < —e, 

T + {lT T * SM int(&))> r +( 7 ^T*SM int (^)) — C ' r — (FT*SM int (&))> T ~( 7 r T*SM int (&)) < ~ € - 

Therefore there exists M such that 

c {(f> s >s) G V Pc : (7 T T *d-SM(.€),V + )g < -e, 

e<s, s < T + (ir T *d-SM(0 ) - e, 

||eii 9 <M}- 

This last set is compact and since Pfi 1 (nfi, 1 (K)) C Pq 1 {ttq 1 {L)) we have that Pff 1 {irf} {K)) 
is a closed subset of a compact set, and therefore is compact itself. Since Pc is a diffeo- 
morphism this implies that ttq is a proper map and thus completes the proof. □ 

Now let x G C°°(SM mt x SM int ). We will eventually need to cut up the operator JVsm 
using such functions. Indeed, by the notation \Nsai we mean the operator whose Schwartz 
kernel is given by x^Ar SM where Kjy SM is the Schwartz kernel of A fsM- With this in mind, 
for the next step we consider compositions of the form 

(x^Sm) 0 7T*. 

As we will see in the next theorem this is still a Fourier integral operator. 

Theorem 2. Define 

L X = (x n sm) o ti** 

Then L x € X-("+D/ 4 ( S M int x M in \ C ' L ; where 

(9) 

C L = {(DF\l£,fj) € T*SM int x T*M int : v e SM int , 3ve SM int 

such that 7 t(v) = it F(v) = F(v), DF\$£ = Dir\lrj}. 

Furthermore, the principal symbol of L x only vanishes at points (DF |* £, fj) for which, when 
v satisfies the requirement in ([9]), x( v ,v) = 0. 

Remark 2. As with Nqm we can easily find an explicit form for the action of L x which is 

L xif}(") = J ( ) x(^>(g)) | d|gj ^-| 1/2 ( 71 /(«)) ds \dSM\ 1/2 . 

Proof. Composing the canonical relation for ir* given by Lemma Q] with that for Msm given 
in Theorem [Tj we obtain Ci x given in Q. If, as in the proof of Theorem [H we show that 
the composition of the canonical relations is clean then this will complete the proof as 
before. 

We refer to points (i)-(iii) given in the proof of Theorem [l] and again show that they 
hold in the present case in which equation © is replaced by 

C = (Cjy SM x C n *) n ( T*SM int x A (T*SM) x T*M int ). 



12 


SEAN HOLMAN, GUNTHER UHLMANN 


First, similar to the previous proof, since Dir |* is injective for all v, and by the hypothesis 
that there are no self-intersecting geodesics, ire is injective and so (iii) holds and the excess 
is e = 0. The statement concerning vanishing of the principal symbol will follow from the 
clean composition calculus as well, and so all that remains is to establish (i) and (ii). 

To begin proving the other parts we introduce a sub-bundle V of T*SM with fibres over 
each point v G SM defined by 

V v = Range(Dvr|*). 

For every v, Dirfa is a linear isomorphism from T*^M to V v , and we denote the inverse by 
(-D7 t|*) -1 : V v —> T*,jM. When these maps for each fibre are combined together the result 
is a bundle mapping which we will write as Dir^ : V —>• T*M. In order to parametrize C 
we introduce an embedded submanifold O of V defined by 

(10) 0={ceV : Dir-^MO) = o}. 

Using the same mapping G as in the proof of Theorem [U let us now show that O = 
Range(G) n V. Indeed, if ( £ V, then 

CM’(MC))) = D^CiDirXZiirviO)) 

= DTT~ t C( 7^(0 (°)) 

= £br”*C(-w(C))- 

Since ( £ Range(G') if and only if {iry (£))) = 0, this proves the claim that O = 

Range(G) n V. 

Now we use the set O to parametrize C via the map 

P C : Vp c = |(C,s) GOxI : T_(7Ty(C)) < S < T+(7iy(C))} -t c 

defined by 

Pel C,») = s)l!f ( „(c ) ,, ) C.C, C, Dn-%) 

where we recall that T is the geodesic flow in SM. Note that since for any v £ SM, and 
any s where the right hand side is defined, 

F(v) = F(V(v,s)), 

we have that 

( 11 ) DFl^DV^syiDFW,^, 

and so _D\h(-, s)|* preserves Range(G). Thus we see that Pc is an injective immersion into 
Range(G) x Range(G) x Range(G) x T*M. To show P^ 1 is continuous we note that it is 
the restriction to C of the map 

pext -1 . j^ an g e ^j x Range(G) x Range(G) x T*M —> Vp c 
P* xt ’~\A,B,C,D) = (B,GJ 1 (B)~GJ 1 (A)). 


defined by 
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Thus P^ 1 is continuous, and so Pc is a smooth embedding and C is a 3n — 1 dimensional 
embedded submanifold. 

To complete the proof of (i) in this case note that Cj\f SM x CV* can be parametrized by 
the mapping 


Pcm sm ^ : Vpc„ sm%c ^ = e T*d.SM xVxRxl : 

0 < s < T+(0, 5 - T+(f) < s < s} 

-A- T*SM int x T*SM int x T*SM int x T*M int 


defined by 


Pc 


A SM 


<C n ,(C,V,s,s) = (DF\ 


T* SM 


«)(5- s ) 


-s&DF | 




C) 




r T*SM®(* 


^11 

77T. 


T* SM 


«)(*) 




^ n T*SM K) 


m (,C,Dn-‘c) 


We now follow precisely the same reasoning as in the proof of Theorem [T| Indeed, if Co € C, 
and X £ {Cj^sm x CV«), then there is a smooth curve a i-A (^(a), C(a), s(a), s(a)) such 
that Pcx SM xC„» (£(0), C(0), s(0), s(0)) = Co and the derivative along this curve at zero is 
X. As in the previous proof, if X £ T C0 (T*SM vnt x A (T*SM mt ) x T*M int ) as well then 
in any local coordinates 


«°) = 5S 


< 2=0 


From this, noting also that 

VsM©( S ' _ S ) ~ ^ i K DF ^n TtSM (i)( S ) 

we conclude that 


«)• 


—s 




<2=0 


Pc^ SM xC n A^(oi),C{a),s(a),s(a)) = — 


Pc(C(a),s{a)) 


< 2=0 


and so X £ T Cq C. This completes the proof of (i) in this case. 

The proof of (ii) follows exactly the same procedure as in Theorem [TJ and we will not 
repeat the details here. □ 


We finally wish to complete the analysis of J\f ( p by considering the composition 

( 12 ) 7 T* O 

We comment that, following from the previous remarks, it is not difficult to show that A 
is given explicitly by the formula 


[ m ( 

( rr + (v) 

JS X M ' 



j TT7?(7^('S)) </>(>(«)) ds | \dS x M{y)\ 1/2 . 


\dv g \ 1 / 2 


The composition of the canonical relations in this case is not as simple as it was in the 
previous theorems. Indeed, the composition in general has multiple connected components 
with one component giving a pseudodifferential operator, and the others occurring only 
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when there are conjugate points. Under some additional hypotheses the extra components 
which appear in the case of conjugate points give rise to Fourier integral operators, whose 
canonical relations and orders can be calculated. This is shown below in Theorem|5]which is 
our main result. For the statement of Theorem[5]we use the following definition, notations, 
and lemmas which also clarify the geometric structure of the operators involved. 

Definition 3. (Conjugate pairs) A pair of unit tangent vectors (v,v) £ SM x SM is 
called a conjugate pair if v and v are both tangent to one geodesic, and 

K (v,v) = ker^ZUr^ o I> w ®| (l , lT+(w) - r+( 0 )) ) p| ker ± {0}. 

If the dimension of Kt v $\ is k > 1, then (v,v) is a conjugate pair of order k. We define 
the set Cr of pairs of regular conjugate vectors to be the subset of SM int x SM mt such 
that given any (v, v) £ Cr there is a neighborhood U of (v, v) in SM mt x SM int such that 
any other pair (w, w) £ U which is also conjugate is conjugate of the same order. Further, 
Cr, k will be the subset of regular pairs of conjugate vectors having order k. Finally, the 
set of pairs of conjugate vectors Cs which are not regular will be called the set of singular 
conjugate vectors. 

Of course it is standard to define conjugate points along geodesics in terms of vanishing 
Jacobi fields, and in fact Definition [3] is equivalent to this. We use Definition [3] because 
it is more convenient in the sequel. In the following lemma we prove that Definition [3] is 
equivalent to the traditional definition of conjugate points. 

Lemma 3. If (v,v) is a conjugate pair of order k (as per Definition \Bf), then there is a k 
dimensional space of Jacobi fields along 7„ vanishing at both 7 t(v) and n(v). Conversely, 
if Jv passses through n(y) and there is a k dimensional space of Jacobi fields along j v 
vanishing at tt(v) and then (v,v) is a conjugate pair of order k. 

Proof. Suppose first that (v, v) is a conjugate pair of order k, and take 

X £ ■,) = ker (jJir\v o P| ker (dtt\ v ^ . 

Suppose that a is a curve in SM such that 

X = d(0), 

and that we have some natural coordinates (x J , v 3 ) on TM such that in these coordinates 

d 


X = a 3 


diJ 


Using X we define a vector field J along 7„ by 

(14) = Dnlyfat) o D v VX. 

Clearly J vanishes at n(v) and n(v), and J is a Jacobi field since it is the variation field of 
the variation through geodesics defined by 

(15) 7(t,s) = vr(T(a(s),t)). 
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Finally, to check that the mapping X H > J is injective (and so there is a A: dimensional 
space of vanishing Jacobi fields), a calculation in coordinates, using the fact that J vanishes 
at n(v), shows that 

(16) V v J(n(v)) = 

Here V is the Levi-Civita connection given by the metric g. Thus the map from X to 
is injective, and so the map X i-a J is also injective. 

Conversely suppose that j v passes through and there is a k dimensional space of 

Jacobi fields along j v vanishing at ir(v) and ir(v). Let J be any such vanishing Jacobi 
field and let 7 (t, s) be a variation through unit speed geodesics with variation field J. 
The geodesics in the variation can be chosen to be unit speed since J must vanish at two 
different points. Then define 

x =£ (o - o)eTsM 

so that in particular X = d(0) where a(s) = ^(0, s). Then (1151) holds, and from this we 
can then see that X and J are related by (1141) . This completes the proof. □ 

We now study the structure of the sets Cr^ by adapting the method of m to prove 
the following theorem. 

Theorem 3. For each k the set Cr^ is an embedded 2n — 1 dimensional submanifold of 
SM int x SM vnt . Furthermore, the set 

J R ,k := {((v,€f), (X,X)) G T(SM int x SM int ) : (v,v) G C*,*, X G K (vfi) , 
forms a smooth vector bundle of dimension k over Cr^- 

Proof. The proof of the first point is an extension of methods found in [30]. Indeed, let 
us begin by considering the exponential map exp(u;) = 7 ^( 1 ) : TM int D TWp M int - 
Let Dpex p be the differential of this map restricted to the tangent space of the fibres. 
Note carefully that we are considering here the exponential map acting on TM mt , but this 
differential is only in the fibre variables. 

Now suppose that (v,v) £ Cn.k- By the results in |30j there exist coordinates on a 
neighborhood W of w = (t + (v) — t + (v))v £ TM int and a neighborhood V of n(v) = exp (w) 
such that the derivative in the radial direction of the A: — 1st elementary symmetric function 
applied to the eigenvalues of Dpex p expressed in these coordinates, which we will label 
<7fc_i, does not vanish. This implies that < 7 ^ 7 ( 0 ) corresponds, through the coordinate map, 
with an embedded 2n — 1 dimensional submanifold of T> eX p \ {0} contained in W with the 
relative topology. Note that ojT^O) is precisely the set of vectors in W with conjugate 
points of order k or higher in V. Now let U be a neighborhood of (v, v) in SM mt x SM int 
such that all other conjugate pairs in U also have order k and shrink U if necessary so that 
for all (w, w) £ U, Xw £ W for some A £ M + . 
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Next consider the smooth map 

p : V exp \ {0} -> SM int x SM mt 

dehned by 



From our construction we have that 

P (a^ 1 (0)np-\u)) = c R , k nu. 

Also, p\ a -i (o)np-!(C/) invertible with inverse defined on its image by 

G(v,v) = (t+(v) - t+(v))v. 

This formula shows that the inverse may be extended to a continuous function on all of 
SM int x SM int and therefore that p\ c - 1 ( 0 )n P - 1 (r/) a t°P 0 l°gi ca l embedding. Since we 

have already shown that aff | (0) np _1 ([/) is an embedded submanifold of T> exp \ {0} we 
also have that p| CT -i (o) n P ~ 1 (U) a smo °th mapping and thus combining all of this we 

conclude that Cp^ n U is an embedded submanifold of SM mt x SM mt . Since ( v , v) £ Cp^ 
was chosen arbitrarily this implies that in fact Cp & is itself an embedded submanifold of 
SM int x SM int . 

For the second point we simply comment that Jp^ can be identified as the sub-bundle 
of another vector bundle Vpj. over Cp t k obtained by pulling back the vertical sub-bundle 
V of T(SM int x SM mt ) ( V = ker(Z?7r) x ker(T>7r)). With this set-up, .Jpj. is the kernel of 
the bundle mapping 

Vp, k 3(X,X)^X- D v * |( W , T+ ( W )_ T+(C ))X G T(SM int x SM int ). 

Since this mapping is constant rank k the kernel is a sub-bundle of Vp k as claimed. This 
completes the proof. □ 

Now we introduce a bundle mapping defined on Jpj. that will be used to describe the 
canonical relation of the normal operator. 

Lemma 4. For each k there is a bundle mapping 

C k ■ Jp,k -»■ T*(M int x M int ) = T*M int x T*M int 

defined by the requirement that for ((v,v), (X,X)) G Jp^ 

( }Di SM \vXUg,i DisM \_xUg ) = D^T(Mint xM int)\\v t v)Ck{((v,v), (X, X))) . 

Proof. First we note that for any v £ TM mt , the kernel of Dtttm\v 5 which is the tangent 
space of the fibre, is a Lagrangian subspace of T V (TM int ) with respect to the symplectic 
form uj g . This can be seen from the fact that \> a is a bundle isomorphism. Indeed 

TTTM = Kt*M ° bg => DtTtm\v = DtXt*M l\, g (v) ° D\>g\ v 

implies that D\> g \ v maps the kernel of Dtttm\v into the kernel of ( v )- Since D\> g \ v 

invertible, and these kernels have the same dimension, in fact D b 9 |„ restricted to the 
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kernel of Dittm\v is an isomorphism onto the kernel of DirT*M\\> g ( v ) which is a Lagrangian 
subspace with respect to the canonical symplectic form u Since \> g is a symplectomorphism 
by definition of u g it follows that }<jsy(D'ktm\v) is a Lagrangian subspace. This point can 
also be shown using coordinates and the formula (flTl) below. 

Now, for ((v,v), (X, X)) £ Jr *, we have that 

Dtttm\v ° DIsm\v X = 0 and Dtttm\v ° DIsm\v X = 0. 


The remainder of the proof is the same for the two cases, and so we just consider the first. 
Indeed, let {x l , v 3 ) be natural local coordinates on TM in a neighborhood of v. Calculation 
in these coordinates reveals 


d 


d 


dvi 


[ x Dis M \vXUg\ ( a% Q x i + ) ~ u g ( Dism\v X, a 1 + V 


d 


dx l 


d 

dvi 


— u g (dism\v X, a Q x i \ ■ 


In the second equality we have used the fact that ker(DiTTM\v) is Lagrangian and Dism\v X £ 
ker(DiTTM\v)- From this we have 


1 DisMlvXUg - u g 


(yDi S M\v X 



dx* £ image(L)7TTM It)- 


This shows that a map is defined by the requirement given in the lemma, and in fact from 
this last formula we can see that it depends linearly on the fibre variables in and 
smoothly on all the variables. Therefore the proof is complete. □ 


It is worth noting that the map Ck can also be expressed in at least two other ways, and in 
fact has a rather simple expression in coordinates. Indeed, if (x l ,v l ) are natural coordinates 
on TM, 

■ d 

Dism\v X = 

and 

H d 

v = 


dvi 


then in coordinates we can calculate 

[dgu 


(17) 


Un = f —— dx 3 A dx l + gn dv J A dx l 
d x j 


From this we find that if ( 77 , 77 ) = Ck((v,v), (X,X)), then in the corresponding coordinates 
x 1 on M 

77 = a l g ik dx k . 

With this in mind, we can see that 77 is also obtained by applying the inverse of the vertical 
lift to DisM\vX, followed by \> g . Another equivalent coordinate invariant definition is to 
take the Jacobi field J defined in (ITT]) . and then 77 = \> g X v J(ir(v)). 

We will prove one more lemma regarding the structure of conjugate points, which is 
really the geometric heart of the proof of our main result. 
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Lemma 5. Let v and v E SM int , r] E T*,sM, fj E T*^M, £ E Tp^d-SM be such that 
F{v) = F{v), 

(18) Dnllr] = and Dir\lr] = DF\l £. 

Then rj(v) = 0, 77 (D) = 0. If v / v then the pair (v,v) is conjugate, and if (v,v) E Cr^ 
then (rj,fj) E Ck(JR,k)- Conversely if ( 77 , 77 ) E Ck(JR t k ) then there is a f E Tf^d-SM such 
that (HID holds. 

Proof. We first observe that since F(j v (s)) is constant with respect to s 

DF\ v (jy v (0)) = 0 

and combining this with the hypotheses we see that 

V(v) = Dn |* 77 ( 7 „( 0 )) = = £{DF\ v (%(0))) = 0, 

and similarly 77 (D) = 0. This proves the first assertion. 

Next assume that v / D. Since F(v) = F{v) if we set s = t + (v) — t + {v), then T(i;, s) = D. 
Therefore 

(19) F(v) = F(^(v, s)) => DF\ V = DF\y o 
and so for £ E Tf^d-SM 

DF\ t v Z = D v *\\ V ' S) oDF\lZ 

which implies 

(20) DttUt] = A^lJ^) o DirWr). 

Now recall that %sm '■ SM mt —> TM int is the inclusion mapping. Then we have 
7 r = tt TM o i S M and s)) = T( ism(v ), s) 

which implies 

(21) Dir}*, = DisM^v ° DtttmII ancl \ v ,s) ° Di SM\l = DisAtll ° D v ^\ 

Therefore, using the fact that the kernel of DirmH is the span of the differential of w i-A 
\w\g, ([ 2 QD implies that 

= F>v^\\ V}S ) ° DtttmIIv + « d 

for some a. Applying the covectors on each side of this equation to a radial vector r based 
at v we obtain 

0 = fj ^Dtt T m\v o + «■ 

Since Dtttm\v ° D v ^\( V:S ) r is parallel to D the first assertion in the lemma implies that 
a = 0 and therefore we conclude that 

( 22 ) DittmIIv = A ^|(„ )S ) ° D-ktmH v- 

We now claim that this implies that (u, D) is a conjugate pair. 




MICROLOCAL ANALYSIS OF THE GEODESIC X-RAY TRANSFORM 


19 


To prove this claim we will set 


We first show that X is in the range of Di$M- To do this, take any Y £ ker(DiiTM\v)- 
Then 

0 = V (Dtttm\vY) = DTTTMllviY) = Ug(X,Y). 

As part of the proof of Lemma |4] we showed that ker(DnTM\v ) is Lagrangian, and so this 
implies that X is in Next take a natural coordinate system ( x l ,v l ) on TM 

in a neighbourhood of v such that g Vj = 5ij at vr(u), and 


(23) 


v = Dtt tm \ 


d 


v dx 1 ' 


Then since X £ ker(_D-7r7\w k>)> in these coordinates 


for some coefficients a?. We thus have 

0 = r)(v) = DtttmIv'T) 


(—) = 


= U) n 



Using dim we finally see that 

0 = a 1 , 

and so X is in the range of D%sm- Thus we can define Di$M Iv 1 ^ £ TSM int . Now we will 
show below that 

(24) X £ ker^Dir T M\v ° Aj^lfos)) f] ker ■ 

By ([2D , after taking transposes, (f2~il) implies that Diqm\v^X £ K v ^ thus proving that 
(v, v) is a conjugate pair. 

Let us begin proving (1M1) . We have already shown above that X £ ker (Dtttm\v)- To 
finish, take any Y £ ker(DirTM\v )• Then similarly to before using also (1221) in the second 
equality we have 


(25) 


0 = DtttmIIv (r) = DittmIvV (A^I^T) = Ug 


X,D„* 


-1 

(v,s) 



Finally, since o 4 $(d, s) is a symplectomorphism we have 


(26) 


o = ujg (a%, s ) X,y) . 


Therefore D v $f\r v>s \X is in the kernel of Dtttm |f which completes the proof of (124(1 . Note 
that we have actually shown that (Dzsm l ^ 1 X, DismIv 1 ° X) £ K( v ^\, and so 

noticing also that 

Dtt T m\Iv = '■xu gi and DutmIIv = t) x u 9’ 

(the second equality follows from a calculation quite similar to ()25|) and (|261l ) we see that 
if (v,v) £ C Rtk , then ( 77 , 77 ) £ C k (J R)k ). 
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Now we assume that ( 77 , 77 ) = C k ((v,v),(X,X)) for some ((v, v), (X, X)) £ To 

complete the proof we must show that there exists £ £ Tf,sd-SM such that (fTKI) is 
satisfied. We split the proof into several steps starting by showing that 77 ( 77 ) = 77 ( 77 ) = 0. 

For this first step we take the same natural coordinate system (x l ,v l ) on TM as described 
above (1231) . Then since X £ ker(Z?7r|„), in these coordinates 

n n 

Di SM \v X = J2 aJ g- 

3 =2 

for some coefficients a?. We have 

V (v) = Dir\lri (J^j = \Di SM \ v xu g (lHsm \v X , . 

Using (11711 this gives proves that 77 ( 77 ) = 0. The same argument shows that 77 ( 77 ) = 0. 

Now, from the first paragraph of this proof we see that (3 £ range(UT|( J ) if and only 
if f3(%(0)) = 0. Since -D7r|*77(7„(0)) = 77 ( 77 ) = 0, we thus conclude that there exists 
£ £ T*^d_SM such that 77 = DF\ f v £. The same reasoning shows that there exists 

£ £ Tp^d-SM such that DnUfj = DF |~£. It remains to show that £ = £. 

For this last step, we begin by noting that a lengthy but routine calculation similar in 
flavor to (|25l) . and making use of the fact that X = s ) X, shows that 

D*U = D v *\\ Vt8) o D7T\lf, =7 DF\l £ = DM (v , s ) 0 DF\l £. 

Finally using (fT9l) (tranposed) this implies 

DF\l £ = DF\l £, 

and so the injectivity of DF |* completes the proof. □ 

We have now reached our main theorem characterizing the structure of the normal operator 
Af. In the statement we use the term local canonical relation to describe the canonical 
relations Ca*. ■ This means that they are each the image of a constant rank map 717 . : C k —> 
CA k such that for any c £ C k there is a neighborhood U of c such that 7 T k (U) is a canonical 
relation. 

Theorem 4. Suppose that Cs = 0- Then the sets 

C Ak = C k (J R , k ) C T*{M int x M int ) 

are either empty or are local canonical relations. On the level of operators, if (p £ C°°(SM) 
is greater than or equal to zero everywhere and Af^ is defined by (ffij) . then we have a 
decomposition 

n—l / M k 

A, i = T + E E a *. 

k =1 \m=l 

where T is a pseudodifferential operator of order —1, and for each k either 
A, (Z j-{n-k+l)/2 (Mint yfint n' . Q 1 / 2 \ 
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where CA km C CA k for each rn, or Mk = 1 and i = 0 if CA k = 0 - Furthermore, T is 
elliptic at every point rj £ T*M int such that there exists a v £ SM mt with rj(v) = 0 and 
cfiv) / 0. 

Remark 3 . In dimension two , there can only he conjugate points of order one and so 
this theorem covers all possibilities in that case. However in dimension three or higher the 
generic case includes singular conjugate pairs. Indeed, according to [I] and m, using the 
notation of [I], singidarities of type D4 occur generically in the exponential map in three 
dimensions or higher, and correspond to singular conjugate pairs of order two. 

Proof. The clean composition calculus does not apply directly to (fT2]l . but as we will see we 
can apply a partition of unity so that the calculus applies to each of the separated pieces. 
Let us begin by looking at the composition of the canonical relations for the operators in 
(fT2l) using Lemma [D and 0 

(27) 

cv* O C L = |( 77 , 77 ) £ r*M int x r*M int : 3 (v,v) £ SM int x SM int 

such that F(v) = F(v), 3 f £ T* F{v) d-SM int 

such that = DF\ f v £, Dn\~ri = DF\~ ^j. 

We will split this composition into different pieces corresponding to different orders of 
conjugate points. First we note that one piece is the diagonal 

A = {(77,77) £ T*M int x T*M int } 

which can be seen to be contained in CV, o Cl by taking v = v such that 77(17) = 0 in ( 1271 ) . 

Now suppose we have (77,77) £ CV* o Cl with 77/77. Suppose v, v and £ satisfy the 
requirements in (1271) . Then since DF |* and Dir]^ are injective we must have v / v and we 
may use the result of Lemma [5] and the hypothesis Cs = 0 to conclude that 

( 71 — 1 

U Ca * 

k=\ 

This union may not be disjoint in general, but we will show that we can still decompose 
the operator Ms into a sum of FIOs each having a canonical relation contained in one of 
the sets in the union. 

We will take the convention that Cr iU = {(f,u) £ SM int x SM mt }. By the hypothesis 
that Cs = 0 and the definition of Cn.h there exists a collection of open subsets Uk of 
SM x SM with disjoint closures such that for each k = 1 to 77 

Crt,k C Uk- 

Thus it is possible to construct a partition of unity {4>k}k=o 011 x SM such that 
supp(0fc) C Uk for k = 1 to n and 

CH- 



supp(0o) C (SM x SM) \ 
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Since this is a partition of unity 

n 

■^<t> ^ 7r * ° ■ 

k =0 

By our construction, for each k the wavefront set of the operator is contained in 

the set 

C' L k = {(DF |* £, -fj) £ T*SM int x : u € 5M mi , 3 5 £ SM rni 

such that (u,u) € supp(^), ir(v) = 7 t T *m(v)i F ( v ) = F (v), DF\l% = DTr\~fj}- 

Therefore, when we microlocalise in the proof of the clean composition calculus (see [10]) 
we obtain smoothing operators except near points in C nf o Ci,,k- For A: = 0 this composition 
is empty, and so we obtain a smoothing operator in that case and this is included in T. 
For k = n the composition is the diagonal A, and finally for k = 1 to n — 1 we have 
C-k, o CL,k = F A k - We will now show that each of these compositions is clean. 

We’ll first consider the composition C nt o Cl , n = A. Based on our constructions so far 
using the results of Lemma [5] we have the following 

C n = (CV, x C Ltn ) p| ( T*M int x A (T*SM int ) x T*M int ) 

= {(77,^1*77,^1*77,77) : v £ SM int , V eT : { v ) M int , r ,( v )= 0 }. 

Using the fact that O introduced in (fTOl) is an embedded submanifold, we easily see that 
C n is also embedded submanifold of dimension 3n — 2 . Furthermore we can see that the 
intersection is clean in the following way. Suppose that <5 = (<5i, 82 , 83 , 64) £ C n and 

{Y 1 ,Y 2 ,Y 3 ,Y\) £ Ts (C^ x Cl,u) P T$ ( T*M int x A (T*SM int ) x T*M int ) =: D n 

where we are considering this set D n as a subset of 

T Sl (T*M int ) x T 52 {T*SM int ) x Ts 3 (T*SM int ) x Ts 4 (T*M int ). 

Then we must have I 2 = Y^. and based on this and an examination of Cl,u we find that 
in fact I 3 and I 4 are determined by Y 2 . Therefore the dimension of D n is at most the 
dimension of CV* which is 3n — 1. However, we can observe that the requirement r)(v) = 0, 
which holds because of the requirement DF \= Dir^f} in the definition of Cl,u an d 
the first part of the proof of Lemma [5j eliminates one more dimension and so in fact the 
dimension of D n is at most 2>n — 2. Since necessarily T§C n C D n and the dimension of 
TsC n is 3n — 2 we conclude that the intersection must be clean with excess e = n — 2. It 
is also easy to see that the projection map ir n : C n -A T*M mt x T*M int is proper and the 
fibres of this map are exactly 

TTn 1 ^) = r|* 77 ; L>vr|*??,??) : r](v) = o} 

which are certainly connected. The clean composition calculus thus implies that ( 7 r* o 
F {<j>®<t>)<t>n a pseudodifferential operator of order (1 — n)/4 — (n + l)/4 + (n — 2)/2 = —1 
which is T, modulo a smoothing operator, in the decomposition given in our theorem. We 
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comment that the excess of the clean intersection corresponds with those v € SM int such 
that rj(v) = 0, and so a principal symbol for T can be found by integrating the product 
of symbols for 7 r* and L^ ( g>^) 0 n over this set. This agrees with formulae which have been 
found for a principal symbol of the normal operator in the past for simpler cases, and 
proves the statement in the theorem on the ellipticity of T. 

Now we turn to the compositions corresponding to conjugate points. If any of the 
compositions C o C kk are empty, then the corresponding operator is smoothing, and so 
can be absorbed into T. If the composition is not empty, then as in the last case we have 
based on Lemma [5] that 

C k : = (C*. x C Lyk ) p| ( T*M int x A (T*SM int ) x T*M int ) 

v), (X, X)), DIsmIv i Di 8 M \vXWg, Di S M ll 'm SM \-xUg, C|((u, v), (X,X))) 

: ((v,v),(X,X))eJ R , k 

Here C k and C k are the first and second components of C k and %sm '■ SM mt —> TM mt is 
the inclusion mapping. This is an embedded 2n + k — 1 dimensional submanifold which is 
parametrized as shown in the previous formula by J Rtk since the mapping 

jR,k 3 ((v,v),(X,X)) (Disult ioisMUxUg, Di SM || ie T*(SM int XSM int ) 

is an embedding. To see that this is an embedding, note that J Rjk is a sub-bundle of V r k 
defined at the end of the proof of Theorem [3l and that the mapping extended by the same 
formula to V R)k is a diffeomorphism onto the bundle O. defined by (fTOl) . pulled back to 
C R k . To see that it is a diffeomorphism, note that it is a bundle map over the identity, 
and then check in coordinates that it is invertible in each fibre. 

We can see that the intersection is clean by a dimension counting argument similar to 
the previous case. As before suppose that 5 = (di, 62 , S 3 ,5/C) € C k and 

(Hi,L2,13,14) e T S (A (7r * )t x A L, k )f]T 5 ( T*M int x A (T*SM int ) x T*M int ) =: D k . 

Just as before we have Y-i determines I 3 and I 4 . In this case we lose n — k dimensions from 
dim(C' Ti ,) since for fixed (v,v) the set 

{DismII ' x D isM \ v xU g : {(v,v), (A, X)) G J R)k } 

is a k dimensional vector space contained in 

{DvU ■■ V G T* <v) M int }. 

Thus, again keeping v fixed, r] is restricted to be in only a k dimensional subspace of the 
n dimensional space T v M mt . Therefore the dimension of D k is at most 3n — 1 — (n — k) = 
2n + k — 1 which is the dimension of T§C k and so as before we conclude that the intersection 
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is clean with excess k — 1. The projections 

v r k : C k -A- T*M int x T*M int 

are proper because the inverse images of sets bounded away from the edge of T*M mt x 
T*M mt are bounded away from the edge in the larger space. This implies that CA k is a 
local canonical relation for each k (provided CA k / 0 )- 

To apply the clean composition calculus as given in [10] it is necessary that the maps 
ir k have connected fibres. However this is not true in general. Nonetheless by what we 
have already done the sets CA k , are local canonical relations. Thus the clean composition 
calculus still may be applied if local representations of the relevant operators are used and 
the decomposition in the statement can be achieved. □ 

We have attempted to make Theorem 0] as general as possible and in so doing sacrificed 
some clarity in the statement. By a slight modification of the proof we could have the 
following simpler corollary which still covers many cases of interest. 

Corollary 1. Suppose in addition to the hypotheses of Theorem [|] that there are only 
conjugate pairs of order 1 in SM mt x SM int , and no two points are conjugate along more 
than one geodesic. Then we have a decomposition 

Jf* = T + Ai 

where T is a pseudodifferential operator of order —1 and 

Hi € T~ n/2 (m** x M^,C A T,n%L xMint ) • 

Proof. The proof is the same as the proof of Theorem 01 except we note that the additional 
hypothesis that no two points are conjugate along more than one geodesic implies that 
the map 7 Ti : Cj -A Ca 1 is injective, and so Cai is a canonical relation. Thus the clean 
composition calculus applies without further decomposition as was necessary in Theorem0J 

□ 


4. Application to inversion of the normal equation 
We now turn to the problem of inverting or given g solving the normal equation 
(28) M^[f\ = g 

for /. To improve the notation going forward let us write 

n -1 / M k \ 

A = Ak -* ) 
fc=l \ m =1 / 

so that with the decomposition from Theorem 0] equation (I28j) becomes 

T [f] + A[f}=g. 

One standard approach is to find an appropriate function space on which this equation is 
of Fredholm type. Indeed, if A is lower order than T in some sense, we may expect this is 
possible. In the setting of Corollary 0| when the dimension n is at least 3, A = A\ will be 
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lower order as an FIO, and if Ca 1 is a local canonical graph then in fact A has appropriate 
mapping properties. It is not always true that Ca 1 is a canonical graph, see examples in 
E3, and a more detailed study of cases in which Ca 1 is not a canonical graph may be an 
interesting direction for future research, but we will say no more about it here. 

To move forward we will need to take ( M,g ) to be a smooth extension with convex 
boundary of (M, g) also satisfying Assumption |T] Precisely this means M d M mt and 
g \m = 9 - Such an extension can always be found, and related to the extension we have the 
restriction maps R : H s (Vi~ int ) —>• H s {^1 ^ int ) which are continuous for s > 0. The adjoints 

of the restriction maps are the extension-by-zero maps i : H s (Q^ irit ) —> H s (Cl~. nt ), which 
are isometric embeddings for s < 0 by duality. In view of the comments in the previous 
paragraph we make the following additional assumption about the extension ( M,g ). 

Assumption 2. Assume the dimension n is at least three, that all conjugate pairs in 
SM int x SM mt are of order 1, and that Ca 1 (see Corollary{J\) is a local canonical graph. 

It should be possible to construct an extension satisfying this assumption when the same is 
true up to the boundary of SM x SM. Thus the assumption is eliminating the possibility 
that there are singular conjugate pairs, or places where Ca 1 is not a graph over dM x dM. 

We will also need to take an intermediate extension M\ satisfying the same requirements 
as M and such that M <£ M\ <s M. We will make use of all of the same objects defined 
on M, also defined in the analogous manner on M and M\, although we will add a tilde 
for objects on M, and also a subscript 1 for objects defined on M\. Thus for example 
F : SM mt -A- d-SM will be the defined in the same way that F was defined just above ([Tj), 
but with M replaced by M. In the same way we have the mapping Tj : SM[ nt —> d-SM\. 

Next we state the main theorem of this section. 


Theorem 5. If Assumption [H is satisfied, cf £ C°°(SM ) is greater than or equal to zero 
everywhere and for every g € there exists a v € SM int with g(v) = 0 and (f{y) 0, 

then the kernel of acting on T 2 (H^ 2 ) is at most finite dimensional and is contained 

in C£°( Furthermore, if F C T 2 (0^ 2 ) is a closed subspace complementary to the 
kernel of then 


(29) 


\X, 






for all f £ T. 


We will provide the proof of this theorem, but first we make a remark concerning its 
significance. The equation (1291) . which shows stability and continuity of from T with the 
iW 1 / 2 norm to T 2 , is intended to match with the hypotheses required for convergence of 
regularizations in Hilbert scales as originally established in m ■ Indeed, using the result of 
m together with Theorem [5] we have the following corollary which shows the convergence 
rate of Tikhonov regularised solutions to the true solution of the problem assuming that 
Xff, is injective. 
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Corollary 2. Suppose that assumption^ is satisfied and X^ is injective on T 2 (fl^). Let 
q > 0 and p > {q — l/2)/2 be given, and suppose /o £ and g £ L 2 (LIq 2 sm ) 

satisfies 

II 9 ~ ^[/olll^^V 2 ^) < e - 

Then for u appropriately chosen the unique solution of the Tikhonov regularised problem 


satisfies 

Wf ~ foWmtitf) ~ 

If the geodesic X-ray transform arises as the linearization of a nonlinear problem, as it might 
in travel time tomography, we also comment that the condition (12911 plays an important 
role in the analysis of convergence of regularisation methods for the nonlinear problem (see 
for example USED). 

We also comment that the use of X^ rather than X allows the result to cover also the 
weighted geodesic ray transform including the limited data case incorporated by setting f 
to be a cut-off function. 

We now turn to the proof of Theorem [5] which will be broken into a series of lemmas. 
The first lemma concerns the continuity of N^. 


Lemma 6. Let Assumption [H be satisfied and suppose that £ C°°(SM). Then for 
-1 < s < 0 

M</> : -> tf S+1 (^M 2 ) 

continuously. 

Proof. Suppose that M\ is an open set such that M ($= M\ (s M int , and take % £ 

such that x( x ) = 1 f° r x € M. Also, let 0 be a smooth extension of <f to M. By Corollary 

|T] using assumption [2] and the continuity properties of FIOs we have 

(30) x"ojvj : -3 H-+\a]L 2 ). 

We also note that from the equation (fT3l) for we can see that for / supported in M 

Rox m o%o X m oHf]=fif 4> [f]. 

Thus we find that for / £ C°°(M) 


mm 


Hs+An 1 ^ 2 ) 


= \\Rox m oNjO X m oi[f]\\ Ha+1{ ^ 2) 
<C\\ X m oN~o X ^ox[f}\\ 1J2 

V M ' 


<C'||i[/]||^ v 2 

V M ’ 

< c||/|| H . (<) 
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where the constant C > 0 may change at each step. The first inequality follows from the 
continuity of R on H s+l { fl~ 2 ) when s > — 1, the second from (1301) . and the third from the 

1 /2 

continuity of i on H S (Q AI ) when s < 0. This completes the proof. □ 

Next we study the continuity of X^ and X^. 

Lemma 7. Make the same assumptions as in Lemma\Q Then 

Xt:H-V 2 $l)£)^L\n y_ 2 5M ), 


and 


continuously 


K ■ L 2 (vTsm) -> 


Proof. Let / £ C°° . Then using Lemma [ 6 ] 


W X 4>if]\\ L2(a^_ 2 SM ) (• A/ ^[/]>/)L2(f2V 2 ) 

^CH/II^^ 2 )- 

As usual, the constant C may change between steps. This proves the first portion of the 
result. The second part follows by duality. Indeed, for / £ C °°(f7^ 2 ) and h £ C°°(Pt-J 2 ^ 

(^[^]) /)z,2(qV 2 ) = (^;^[/]) L 2(nV 2 s 

< II/ZII-,1/2 \11L7J||r2/o^/ 2 


d-SM' 


< c\\h\\ L2 ^ 


d_SMl 

mm 


d_SM 


L 2 Rlf_ SM ) 

H- 1 / 2 (M)‘ 


Taking the supremum over / such that ||/||#-i/ 2 ( M ) = 1 we obtain the result. □ 

We now begin the proof of Theorem [5l 


Proof of Theorem 0 Let d$M = {x £ M : dist(x,<9M) < d}. Then by the result of [29], 
since dM is strictly convex and compact, there exists a <5 > 0 such that if f £ Ker(A^) 
acting on L 2 (fl^ 2 ), then f\g s M = 0- Now introduce open sets Mi, M 2 , and M 3 such that 

d s M c C Mi m M 2 <e M 3 <e M int . 

Also, let xi e C“(M 2 ) and X 2 £ Cf°(M :i ) be such that Xi(x) = 1 for x £ Mi, and 
X 2 (x) = 1 for x £ M 2 . 

By Corollary [T] and Assumption [ 2 ] we have a decomposition 

JVcj, = T + A 

where T is an elliptic pseudodifferential operator of order —1 and A is an FIO of order 
— n/2 whose canonical relation is a graph. 
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Since T is elliptic it is possible to construct a parametrix T 1 which is a properly 
supported pseudodifferential operator of order 1. This means in particular that 


Y -1 o T = Id + 7£ 

where 7 Z is an operator with smooth Schwartz kernel. Furthermore, we can construct Y ” 1 
such that its Schwartz kernel is supported sufficiently close to the diagonal so that 

Xi* o Y -1 = xT o Y _1 o (x ™) 2 

By the hypotheses, constructions above, and mapping properties of FIOs 
X>Ao X ? : iT(< 2 ) -A tf s+3 / 2 (< 2 ) and ■ YY s + 3 / 2 (0^ 2 ) -A H s+1 / 2 (< 2 ) 

are compact for any s. Putting these together and using Sobolev embedding we find that 
X? ° T - 1 o ( X -) 2 oAo X ?: H*^ 2 ) -A H S (Q]/ 2 ) and X ? o R : H\Q$) -A H S (Q]/ 2 ) 
are compact for any s. 

1 /2 

Now suppose that / € H s (Qj^ ) has support contained in M\. Then using the construc¬ 
tion from the previous paragraph we find that 

(31) / = X ? o Y- 1 o ( X -) 2 o Af4f\ - xr o Y" 1 o (x^) 2 o A o xH/] - xT ° «[/]• 

1 /2 1/2 

By the previous paragraph the operator K : H s (Q^ f ) —y H ) defined by 

K = X ?oT- 1 o( X ™) 2 oAo X ™- X ™oR 
is compact. Using also the mapping property of x™ 0 Y _1 o (x ™) 2 we find that 


(32) 


HA<) ~ ° 


'H‘+Hntf) + [/]ll iyS ( jr2 V 2 ) > 

This holds for / with support in M\. For convenience we now introduce the notation 

mn)£) = {/ e H*^ 2 ) : supp(/) C Mi} 

1 /2 

which is a closed subspace of H s (). From (f32l) we can establish using standard methods 
(see for example [25]) that 

ker^n 

1 /2 

is hnite dimensional. Also we can show that for any closed subspace T\ of ) on 

which Ah is injective 

( 33 ) - c mm H s + 1 ( n]i 2 ) 


for all / 6 J 1 and any s <E M. 

Now suppose that / £ Ker(A^) acting on L 2 (fi]^ 2 ). Then since f\e s M = 0 we have that 
the support of / is in M\. On the other hand from (1311) we have 
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from which we may conclude that / £ Repeating this we obtain by a boot 

1 /2 

strapping argument that / £ This proves the first statement of the theorem, 

and it only remains to prove (129(1 . 

Actually, half of (12911 has already been proven in Lemma [?1 and all that remains is to 
prove the stability estimate 

( 34 ) 11^11/1-1/2(0^2) < C\\X,p[f] lli2(nV 2 gM ) 

for all / £ T. For this we apply (1331) with M and M\ replaced by the extensions M and 

- -—- _1 /o i /o -—- 

M\. Then T\ = ifJ 7 ] is a closed subspace of H 0 ' (9~ ) on which X^ is injective since the 

extension-by-zero operator i : iL _1 / 2 (9^ 2 ) -A- iL _1 / 2 (9~ 2 ) is an isometric embedding for 
s < 0. Using this embedding property again we find that 

II^Ih-i/rO = ll i [/]ll J H-v 2 (n_ t/2 ) < c, ll-^ oi [/]ll//i/2 ( oR2 ) . 

Applying Lemma |7] on M to this last estimate we have 

h- 1 / 2 ^ 2 ) - c \\Xj > °'Af]\\ L 2(n 1/2 _)■ 

v M ’ v a_SM J 


For the final step we use the fact that R o J\f^ o i[/] = A//,[/] for / supported in M which 
we have also used in the proof of Lemma [ 6 ] and follows from equation (1131) . Using also the 
fact that R and i are adjoints we have 


IT- 


v d_SM / 


= ^oi\I], Xi oi\I]) LHnW 


d_SM 


= <P 0 ^ 0 i[/]} L 2 (< 2 sM) 
= (f^4>[f]) L 2 { ay* sM) 

= (^■[/]>‘ ; ^[/])L 2 ( n i/ 2 sM ) 

HIWlll* 1/2 r 


Thus (1341) is proven which completes the proof of Theorem [H 


□ 


5. Conclusion 

While we have gone some way towards completing the microlocal analysis of the geodesic 
X-ray transform for nontrapping manifolds, a number of questions remain. Microlocal 
analysis of the two dimensional case is complete, and indeed in EE! it is shown, in the two 
dimensional case, that when there are conjugate points X can actually cancel singularities. 
That is, there exist non-smooth distributions / such that X[f] is smooth. This has the 
consequence that stable inversion is not possible between any Sobolev spaces, and so the 
inverse problem of recovering / from X[f] in this case is severely ill-posed. We have proven 
that this does not occur in three dimensions or higher provided that Assumptions [T] and [2] 
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are satisfied, and indeed that the problem is only mildly ill-posed in that case. When the 
assumptions are not satisfied more work is required to determine the degree of ill-posedness. 

The failure of Assumption [2] may occur in at least two ways. Firstly, as pointed out in 
Remark [3j in dimension three and higher the hypothesis of Theorem [4] that Cs = 0 fails 
for generic metrics, although at least in three dimensions this only happens at isolated 
points which have D 4 type singularities. Study of the normal operator M ( p near such points 
is therefore required for a full understanding of the microlocal properties of X<p in three 
dimenions (which is likely the most interesting case for any application), and in particular 
understanding of whether the inversion is mildly or severely ill-posed. 

Another way the Assumption [2] may fail is the additional requirement that Cai be a 
canonical graph. It is at the moment unclear whether this is satisfied generically, if not 
whether it fails only at isolated points, and what precise impact it might have on the proof 
of stability estimates as shown in section HI The canonical graph assumption is required 
for the continuity of A\ between appropriate Sobolev spaces, but weaker versions of such 
continuity may still hold even when Ca\ is not a canonical graph. 

As mentioned in the introduction, the geodesic X-ray transform for tensor fields is also 
of interest and in fact arises naturally in travel time tomography. It is likely the method 
used in this paper could be extended to the tensor field case with some adjustments, and 
this is reserved for future work. 
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